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' Abstract 

We discuss the processes 77 tt^tt^ and r] — > 7r''77 O(p^) in the mo- 
mentum expansion. The calculation involves tree-level, one-loop and two-loop 
diagrams of a chiral effective lagrangian which is obtained by a bosonization of 
the NJL model. The importance of integrating out meson resonances (reduc- 
tion) is pointed out. Our final results for the total cross section of 77 tt^tt^ 
are in good agreement with the experimental data of the Crystal Ball Collab- 
oration. For the width of the t] — > -^"''77 decay we obtain the value 0.11 eV 
which has to be compared with the experimental value of (0.84 it 0.18) eV. 
Alternatively, taking empirical parameters from a vector-meson-dominance 
model the prediction for the decay width is 0.35 eV. We present a prediction 
for the differential decay probability as a function of m^^/m^. 



* Supported by Deutsche Forschungsgemeinschaft 



The theoretical interest in the reaction 77 tt^-k^ dates back to the seventies 
when predictions for the electromagnetic polarizabilities of the charged as well as 
the neutral pion were obtained in the framework of current-algebra techniques 
and chiral quantum field theory 0. These polarizabilities are a signature of the 
underlying structure of particles, similar to the electromagnetic root-mean-square 
radius, and a large number of different predictions for these parameters has been 
obtained in various models (for an overview see, e.g., Refs. [^). The possibility of 
investigating the 77 tt^tt^ amplitude via the e"'"e~-annihilation process as well as 
the photoproduction in the Coulomb field of a nucleus was addressed in Refs. 0]. 

In the meantime, 77 —>■ tt^h^ cross section data from threshold up to the p- 
resonance region were provided by the Crystal Ball Collaboration ||^. On the theo- 
retical side the framework of Chiral Perturbation Theory (ChPT) ||^, |^ provides an 
ideal tool to systematically study low-energy amplitudes involving Goldstone bosons 
and their interactions with external fields, such as the electromagnetic field. In Refs. 
Ig] the amphtude for 77 vr^vr^ was calculated to O(p^) in ChPT, and the result 
was found to be given entirely in terms of one-loop diagrams involving vertices of 
0{p^). In other words, there are no tree-level diagrams at 0{p^) and O(p^) and 
thus the one-loop diagrams are finite. However, the one-loop calculation in ChPT 
disagrees with the data even near threshold. The inclusion of a Born contribution at 
0(j9^), obtained either from quark loops or from vector-meson dominance, results in 
too small a contribution to yield agreement with experiment 0. On the other hand, 
the application of dispersive methods leads to a considerable improvement since they 
take account of important unitarity corrections corresponding to rescattering effects 
of higher order 0. A full two-loop calculation at 0(/) within SU{2) x SU{2) ChPT 
was carried out in Ref. [|l^. The O(p^) counterterm contributions were estimated 
with resonance saturation and the total result was found to be in good agreement 
up to an invariant mass of 700 MeV. Finally, 77 ix^-n^ was also considered in 
the context of Generalized ChPT up to one-loop order corresponding to 0{p^) in 



this counting scheme |]12 



In the framework of chiral S'f/(3) x S'f/(3) symmetry the decay process 77 7r°77 
is closely related to 77 — -k^-k^. At O(p^) in ChPT the prediction for the decay 
width [13 1 was found to be two orders of magnitude smaller than the measured value 



14|| . The pion loops are small due to approximate G-parity invariance whereas the 



kaon loops are suppressed by the large kaon mass in the propagator. A consider- 
able enhancement was obtained with resonance saturation for some counterterms of 
higher orders in the momentum expansion. In Ref. [|l^ symmetry-breaking terms 



proportional to the quark masses were not considered at 0{p^). Such counterterms 



were, however, included in Refs. [15, M. In Ref. [pH] they were estimated in the 



framework of an extended NJL model [|T^ whereas in Ref. |jT6| the experimental 



decay width was used to fit one of the corresponding coefficients. Finally, a comple- 
mentary approach was used in Ref. where the rj 71^77 decay was calculated in 



a phenomenological quark model using the quark-box diagram. A good agreement 
with the experimental value for the decay width was obtained with a constituent 
quark mass of 300 MeV. 

It is the purpose of this work to present the results of a consistent calculation of 
the processes 77 — * vr^vr'' and r] —>■ 71^77 at 0{p^) in the momentum expansion. Ac- 
cording to Weinberg's power counting scheme the calculation involves tree-level, 
one- and two-loop diagrams. The effective action up to 0{p^) in terms of collective 



meson degrees of freedom is obtained by bosonization |T9[ of the NJL model ^0 
This effective action, in addition to the pseudoscalar mesons, still contains scalar, 
vector and axial-vector degrees of freedom. In order to determine the structure co- 
efficients of the effective chiral lagrangian at 0{p'^) [0 and 0{p^) one has to 
integrate out the meson resonances. The method of superpropagator regularization 
was used in order to fix the UV divergences which for the first time show up at 
0(/). 

We start from the generating functional 

Z = Jv^V^^VVVAexp[iS{^,^\V,A)], (1) 

corresponding to the following action for scalar (S*), pseudoscalar (P), vector (V^) 
and axial-vector (A^) collective meson fields, 

^\V,A) = I d'x [ - ^tr($t$) - -^triV.V' + A.A'^) 



4Gi 
+ log(det(2D)) 



(2) 



This action is obtained by first bosonizing the effective action of the NJL model 
and then integrating over the quark degrees of freedom. In Eq. Gi and G2 are 
parameters which are fitted to empirical input (see Eqs. (|ll|) and (|13|) below for 
details), $ = 5 + iP, and D refers to the Dirac operator 

= + 4n) - + mo)]PR + + 4l) - + mo)^]PL, (3) 

where mo is the current quark mass matrix, Pr/l = |(1±75) are chiral projectors and 
A^'^^ = Vu^Au. The electromagnetic interaction can be included by the replacement 



Vfj, ^ + ieA^Q, where Q is the quark charge matrix, Q = diag(2/3, —1/3, —1/3). 
We express $ using a nonhnear reahzation of chiral symmetry. 



where 



exp 



.V2Fo 



V2' 



TT 



V 



TT 

K- 



(4) 



'Ivs + j=,Vo ) 

represents the pseudoscalar degrees of freedom. Fq is the bare vr decay constant. 
The 3x3 matrix contains the scalar fields and is expanded around its vacuum 
expectation value /x, 

= /i + a{x). (5) 

The constituent quark mass /i is the solution of the gap equation. 

For the processes under consideration, up to and including O(p^), only the even- 
intrinsic-parity sector of the chiral lagrangian is required . This sector is obtained 
from the modulus of the logarithm of the quark determinant and can be calculated us- 



ing the heat-kernel technique with proper-time regularization |2^, ^ . This method 
has been used in Ref. to obtain a prediction for the structure coefficients of the 
general effective lagrangian of O(p^) and 0(p®), respectively 0, The result of 



Ref. 1^ explicitly contains, apart from the pseudoscalar Goldstone bosons, scalar, 
vector and axial-vector resonances as dynamical degrees of freedom. However, in or- 
der to avoid double counting when calculating processes involving Goldstone bosons 
and photons, one has to integrate out (reduce) these resonances in the generating 
functional of Eq. (|I|) and thus one effectively takes resonance-exchange contribu- 
tions into account. As a consequence of this procedure the structure coefficients of 



27]. 



pseudoscalar low-energy interactions will be strongly modified [0, | 

In order to perform the integration over the scalar, vector and axial-vector fields 
in Eq. (|T]) we made use of the fact that the modulus of the quark determinant in 
Eq. (^) is invariant under local chiral transformations of the fields p7| , |28| . This 
allows us, with a specific choice for the chiral transformation (unitary gauge), to 
eliminate the pseudoscalar fields from the Dirac operator, Eq. (^). At the same 
time, introducing $' = $ — itlq and renaming $' — > $ generates the mass term 
for the pseudoscalars from the Gaussian part of Eq. (0). Furthermore, interactions 



between the pseudoscalar degrees of freedom and the transformed vector and axial- 
vector fields are generated in the Gaussian part. The masses of the scalar, vector and 
axial-vector mesons are sufficiently large in comparison with the Goldstone boson 
masses, and thus it is possible to integrate out the meson resonances using their 
respective equations of motion in the static limit. These equations result from a 
variation of the effective action of Eq. (0) by neglecting terms of O(p^) and higher 
in the logarithm of the quark determinant. The remaining part of the action then is 
quadratic in the resonances, in particular, there are no terms containing field strength 
tensors. 

The invariant amplitude = ieie2T^u of the process 7(q'i)7(q'2) ct(Pi)^(P2) 
can be expressed in terms of two functions A and B as 



A{s, V 
+B(s,i^ 



:9fii' — Q2fj.qii 



/ 2 / 2 2\2\ 



+2((z/ + ml- ml)q2^A^ - {u + ml - m^)A^gij, 



(6) 



where s = {qi + ^2)^, = 2pi ■ (g2 - qi), and = (pi - p2)^. The amplitude 
for the process a{pi) — > &(p2)7(q'i)7(q'2) can be obtained from Eq. (|]) using crossing 
symmetry, namely, by performing the replacement qi —>■ —qi andpi —pi- However, 
for the decay channel ri{k) — > 'n'^ {p)'y{qi)'j{q2) , it turns out to be more convenient to 
use the parameterization 



A{xi,X2)[g^,u{qi ■ ^2) - qiuq2fi 



V 



(7) 



where Xi = {k ■ qi)/m'^. 

The prediction for the amplitudes of Eqs. @ and (0) will involve the structure 
coefficients Li of the Gasser-Leutwyler lagrangian in one-loop diagrams at 0{p^) as 
well as new coefficients di from Born diagrams at 0{p^). It is straightforward to 
obtain the effective lagrangian at 0{p^) contributing to the processes under consid- 
eration from the most general representation of Ref . |^ , ^ 



+4^M-^^^tr (q2) tr [d'^Uod^U^,] + deJ^^^mi [Q^) tr {d^Uod'^U, 



^Note that there are different conventions for the definition of the coefficients dj. 



+d,T^^T^hi (d'^U^UlQ) tr (dpU^UlQ 



In Eq. (H), Tp,v = d^j^Au — d^A^ is the ordinary electromagnetic field strength tensor, 

TT (-V^^O. 



= 1 71 + 7i + 71*. "7! + 75 + 7S"»- "l/s"' + 75"° 

and X = diag(xM, xl) Xs) = — 2mo <gg> Fq'^ is the mass matrix, where <qq> is the 
quark condensate. Previous calculations considered the counterterms of Eq. (§) with 
various degrees of approximation. In Ref. [|l^ only single-trace terms in the chiral 
limit were taken into account. In Refs. []TT], |T3| the chiral symmetry breaking term 
proportional to was included and Ref. [1^ also took d^ into account. The double- 
trace terms proportional to d^ - dg typically do not appear in effective lagrangians 
derived from the bosonization of NJL type quark models. 

In the NJL model only the structure constants di, ^2, ^3 contribute to the Born 
amplitudes of the processes 77 tt^tt^ and rj — > 7r°77 at 0{p^), respectively. 



g4g2 [55 

—dis + -d2{s - 2ml) + d^i^xl + Xl) 



9^ 



di 



(9) 



and 



gg2 f g r 



Xu Xd 



6(m2 — 

-\{xl-xl)Qi 



16e2 



Co 



2 + 



(rfi + U2)m\{xx + X2) - 4^2^^ + 4c/3(4x^ + Xd) 

((ii + 4(i2)m^(xi + X2) 
- ^ [d2ml - d^{Axl + Xd + 4x')) } , 

5 xi 



X^ 

3 — 



+ o(x^-x^)ei 



(10) 



In Eqs. ( pUp = cos^^ — v^sin6', where Q = —19° is the 77 — 77' mixing angle. 



rjs = T] cos 6* + 77' sin 6* , rjo = —rj sin + rj' cos 6* , 



and furthermore we have introduced 

(cos^ - V2sm9)^ (sine + ^2 cosy 
m„ — mi m , — mi 

Note that the t] decay amphtudes of Eqs. (|10|) also include contributions of the pole 
diagrams with vr^— 77 and 7r°— 77' transitions. 

We now turn to the determination of the structure coefficients within the frame- 
work of the NJL model. It is a well-known fact that the elimination of the resonance 
degrees of freedom gives rise to a substantial modification of the structure constants. 
At O(p^) such a reduction leads to a redefinition of the decay constant Fq and the 
mass matrix x- To be specific, the identification of the decay constant before and 
after reduction is given by 

^0=^^ — ^o=Za^^, (11) 
respectively, and similarly ^ for x 

X=-2mJl-^^e-^'"^') X = ^, (12) 

where y = r(0,/x^/A^), fi is the average constituent quark mass, A is the intrinsic 
cutoff parameter, and 

zf-iu^r^. nu^- (13) 



AG. 



2 



The incomplete gamma function is defined as r(n, x) = J^dte ^t"" ^. In Eq. 
we have introduced 



4 



*(2.-i; 



487r 



, (m^) =m^(l+7), 7 = ^g^2 



The parameter Z\ of Eq. (|T3|) corresponds to the n — Ai mixing factor and has the 
phenomenological value 



Zl = f ^ [ ^ 0.62, 



1-7 

where we used the following empirical input, rrip = 770 MeV, = 1260 MeV, and 

9v = Qp-K-K = 6.3. On the other hand, with the special choice Z\ = 1/2, Eqs. (^3|) 

and (|Tl|) reproduce the well-known Kawarabayashi-Suzuki relation, = 2gyFQ. 

^Using the gap equation it can be shown that both expressions for x Eqs. ( p^ are equivalent 
for /x^/A^ < 1. 



A full calculation of the it and K decay constants at O(p^) allows to fix the 
parameters y and x = — /iFq /(2 <qq>) for given values of Z\ and ^, by identifying the 
decay constants with their empirical values. In the following we will use Z\ = 0.62 
and fi = 265 MeV, from which we obtain y = 2.4 and x = 0.10. These values 
correspond to Fq = 90 MeV and <gg>i/3= _220 MeV. 

At O(p^) the reduction of the resonances ||2^ leads to the following modification 
of the structure coefficients of the lagrangian introduced by Gasser and Leutwyler 



{Li — jo^h), 



'■4 

k = 2(i 
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h = x{y - 1) 
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^10 = ~77 

6 



ired ~ 

*1 ~ O/l 



1 

24 



1 (v 



ired 

ired 
'3 

ired 

ired 

ired 
'6 

ired 

h 

ired 
(■a 



ired 



ired 

no 



1 

~6 
0; 



Zl + 3(Z: 



^1 



1.08 /i; 



1.54/, 



0; 
0; 



0.60 /, 



5 1 



16 ^ 
1 



1.07/s 



^1 - 



6 



1.12/9; 
1.86/10, 



(14) 



In order to obtain the expressions for the reduced coefficients of Eq. (0), the static 
equations of motion of the scalar, vector and axial-vector resonances have been 
applied. In such an approach scalar resonances can only modify and /§. Note that 
the above results are in agreement with those obtained in Ref. |jl^ except for /g*^'^ and 
/g*^"^ (see Sect 5.5 of Ref. |l^). The disagreement originates in a different procedure 
of integrating out the scalar resonances. We will come back to this point below when 
discussing higher-order corrections to the static equations of motion. 

We will now discuss those structure constants di at O(p^) which do not vanish in 
the NJL model. Before reduction we obtain 



d3 



F2 1 

-^0 



167r2 ^2 24 



Nc 1 



X 



-9.13 X 10" 



1.83 X 10" 



do 



F2 1 
-^0 ^ 



4.57 X 10" 



167r2 12 

The first two constants coincide with the results of Ref. 



167r2 48 

(15) 

. The reduction of meson 



resonances in the framework of applying the static equations of motion generates the 
following modifications 



, ^ _J:!r^EIJ_z\ = -3.51 X 10-^ , dl^' = J%^-Z\ = 1.76 X 10-^ 



AT p2 1 

^ 167r2 /i2 96 ^ ^ ^ 

In this context we note that the modification of the first two structure coefficients 
results from the application of the equation of motion to vector and axial-vector 
resonances. This change amounts to a multiplication of the original coefficients di 
and d2 of Eq. (|T3p by a factor Z\. The situation for d^ is qualitatively different. In 
this case the application of the equation of motion to the scalar resonances modifies 



this coefficient. Let us compare our results for d^ with those of Ref. ||T5|. We 
agree for the coefficients d'['^'^ and d'^'^'^ but differ with respect to rfg^'^. In order 
to understand this discrepancy we note that two different techniques were used to 
eliminate the resonances. In the treatment of scalar resonances the method of Ref. 



T5[| involves operators with derivatives which are beyond the scope of our treatment 
using the static equation of motion. A comparison with Eqs. (23), (32) and (38) of 
Ref. shows that such operators are the origin for the difference in ^2'^'^. However, 
there is another interesting observation. Even though our final expression for rfg*^"^ is 



the same as Eq. (40) of Ref. |]T5l our result originates entirely from the reduction of 
scalar resonances whereas in Ref. |]T3| it is the sum of a scalar resonance contribution 
(see Eq. (39)) and a quark-loop contribution (see Eq. (23)) for which we have no 
analogue. 



Finally, we have also investigated in our approach those results of Ref. |15[ which 
correspond to the inclusion of operators containing derivatives when integrating out 
the scalar resonance. To this end, after a unitary gauge transformation of the modu- 
lus of the quark determinant, one has to keep also higher-order terms in the effective 
action of Eq. (§) which are linear in the scalar field a{x) and which contain the 
coupling to vector, axial-vector fields and field strength tensors. Such higher-order 
terms lead to a modification of the static equation of motion for the scalar resonances 
and thus give an additional contribution to the structure coefficients l'^^'^ and d2^'^, 

4 y 

dr^'--^ = = 1-02 X 10-^ , (18) 



Table 1. Modification of tlie coefficients ai, 02 and b of Eq. ([TP|) due to tlie reduction 
of meson resonances. Af = Nc{4:7!-Fq/ fi)"^ = 54.6, Z\ = 0.62. 



Coeff. 


Without 
reduction 


Reduction of resonances 


Vp- and yl^-fields 
in static approx. 


(T-fiold 


Sum 


Static approx. 


Higher-order correct. 


ai 


'§{Vix~l)K = 12.1 


-izlAr=-15.6 


f,ZlN=12.5 


-izi(i-i)Ar = -9.o 


-12.1 


0.2 


l,N = lb.2 


^Z\M = 5.S 







10.3 


b 




T§8 ^1^ = 0-97 








0.97 



wliicli agree witli Eq. (155) of Ref. |jT3 and Eq. (38) of Ref. [|T5|, respectively. The 
total result for the coefficients l^^*^ and rfg^'^ after reduction of the vector, axial-vector 
and scalar degrees of freedom then is the sum of the contributions of Eqs. (|I4D and 
(|l^) and (|T6|) and ([I8D , respectively. It is worth noting that the considered higher- 
order terms also modify the static equation of motion of axial-vector resonances. 
However, this modification does not lead to any new contributions for either the 
structure coefficients Li or di. 

For the purpose of comparing our numerical results for 77 — -k^tx^ with those of 



Refs. Ipn] , pTSjl , it is convenient to introduce the following parameterization |Tl| of the 
Born contribution at 0{p^) for the amplitudes A and B of Eq. (^, 

A, = + B, = - (19) 

' (167r2F2)2 ' (167r2F2)2- ^'^^ 

The coefficients Oi, 02 and b are related to di, d2 and ^3 by 

= (47r)^H32(d3 - d^) , a, = {AT:Y^-^2{d, + M,) , b = -{A^xY^^d, . 
y y y 

Our results for and b are summarized in Table 1. Clearly, the reduction of the 

resonances leads to a large modification of the coefficients. However, one has to keep 

in mind that the effective action after the reduction describes the interaction of only 

pseudoscalars and photons. Thus the modified coefficients should not be treated 

as additional corrections to the nonreduced coefficients of Eq. (|T3p. A summation 

of quark-loop contributions and resonance-exchange contributions to the structure 



coefficients as in Table 1 of Ref. , in our opinion, leads to double counting. 

Before comparing our values of the 0{p^) structure coefficients with those of Ref. 



Tl[] we provide a prescription for relating results in different renormalization schemes. 



In our approach UV divergences, resulting from meson loops at O(p^), were separated 
using the superpropagator regularization method which is particularly well- 
suited for the treatment of loops in nonlinear chiral theories. The result is equivalent 



to the dimensional regularization technique used in Ref. |Tl]], the difference being 
that the scale parameter /i is no longer arbitrary but fixed by the inherent scale of 
the chiral theory, namely, fl = An Fq. In order to compare the two methods the UV 
divergences have to be replaced by a finite term using the substitution 



(C-l/e) 



C, 



SP 



2C+H-- 



d 



-(logr-(2-- + 2: 



z=0 



where C = 0.577 is Euler's constant, e = (4 — D)/2, and (3 is an arbitrary constant 
resulting from the Sommerfeld- Watson integral representation of the superpropaga- 
tor. The splitting of the decay constants and is used at O(p^) to fix Csp ~ 3.0. 



For our numerical comparison with the two-loop calculation of Ref. |]TT| we made 
use of the parameters Lj and di corresponding to Tables 1 and 2 of Ref. [|lT|. In 



particular, from the numerical values of the parameters Oi, 02 and b of Table 2 of 
Ref. Q 



BGS 



-39.0 



12.5, 



lBGS 



3.0 



one obtains 



df 



BGS 



-10.8 X 10" 



d 



BGS 



4.29 X 10" 



-0.10 X 10" 



(20) 



Our predictions at O(p^) and 0{p^) for the 77 tt^tt^ cross section near threshold 
are shown in Fig. 1. The calculation at 0{p^) contains Born, one-loop and two- 
loop diagrams. In our two-loop calculation only diagrams which are factorizable 
and which can be calculated analytically were taken into account. Two-loop box 
diagrams and acnode graphs cannot be calculated analytically but the numerical 



estimates of Ref. |]TT| indicate that their contributions are small. As was already 
discussed in Ref. [0, the predictions of the NJL model for the coefficients d^'^'^ 



and ^2^'^ are about a factor one half smaller in comparison with the vector-meson- 
dominance model (VMD) (see, Refs. |P, |l5l). The coefficients di and d2 in the VMD 



model can be obtained from Eq. (K^ by the replacement 



^1 = ^ 



6 flGTchvfJ'V' 



Nc \ my 



0.82, 



(21) 



with my = m 
the decays V 



and where the coupling constant hy = 3.7 x 10 ^ is extracted from 



71'j. This has to be compared with the prediction of the NJL model, 
■ 0.62. We have taken account of this uncertainty by 
showing the results for both Z"^ = 0.62 and Z"^ = 0.91. The results of our calculations 



2.5 X 10-2 for Zl 



with the parameters of Ref. |11| are also shown in Fig. 1. Numerically they are in 
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Fig. 1. Cross section for 77 tt^tt^ as a function of the invariant mass W = m^o^o for 
W < 0.7 GeV and | cos^*| < 0.8 where 9* is the angle between the beam axis and one 
of the vr'^ in the 77 center-of-mass system (c.m.s). The data are from the Crystal Ball 
experiment The dotted line represents the one-loop calculation at O(p^). The dashed 
line corresponds to the calculation at 0{p^) without reduction of the resonance degrees 
of freedom. The dash-dotted lines corresponding to two different values of the parameter 
Z\ are a measure for the uncertainty in the reduction of the meson resonances. This 
uncertainty is due to the difference between the NJL prediction and the empirical value 
for the coupling constant hy- The solid line corresponds to the values of the coefficients 



Lj and dj used in Ref. I 111. 



Table 2. Contribution of various diagrams to the rj — > tt 77 decay width. 

C_!^o^= (0.84 ±0.18) eV. 



Amphtudes 


Without 
reduction (eV) 


With reduction (eV) 


Zl = 0.62 


Zl = 0.91 


1-loop 

0{p') 


TTTT-loOpS 


1.3 ■ 10-3 


1.3 ■ 10-3 


1.3 ■ 10-3 


KK-loops 


6.2 ■ 10^3 


6.2 • 10-3 


6.2 • 10-3 


Born 0{p^) 


0.22 


0.11 


0.45 


1-loop 

0(/) 


TTTT-loOpS 


1.9 ■ 10-^ 


6.9 ■ 10-5 


8.6-10-^ 


i^K-loops 


4.1 ■ 10-2 


1.9 ■ 10-3 


2.7-10-2 


2-loop 
0(/) 


TTTT-loopS 


3.2 ■ 10-^ 


3.2 ■ 10-^ 


3.2 • 10-^ 


Tri^-loops 


3.1 ■ 10-3 


3.1 ■ 10-3 


3.1 • 10-3 


KK-loops 


1.4 ■ 10-5 


1.4- 10-5 


1.4- 10-5 


Total 


0.14 


0.11 


0.35 



a good agreement with Ref. |]Tl|; even for 7717^77 as large as 700 MeV the difference is 
only about 7%. 

For the decay width of rj —>■ 7r°77 we obtain after the reduction 0.11 eV and 
0.35 eV corresponding to Z\ = 0.62 and Z"^ = 0.91, respectively. On the other 
hand, using the parameters of Eq. (EOf) one finds 0.18 eV. These results have to be 



compared with the experimental value (0.84 ± 0.18) eV [T^. The contributions of 
different diagrams to the decay width are shown in Table 2. These results clearly show 
the dominance of the Born contribution. It is a well-known fact that calculations 
of the decay width at 0{p^) tend to come out too small in comparison with the 
experimental value [0, |T5|, |1^ . This failure indicates that either higher-order terms 
are required or higher-order resonances have to be included or both. 

Finally, we have also tried to fit the coefficients di, 62 and c/3. However, due to 
a strong correlation between the coefficients di and ^3 it was impossible to find a 
stable minimum from a fit to the 77 vr^vr^ cross section and the 77 71^77 decay 
width. The strong correlation is related to the fact that the m-,;-,^ dependence of the 




Fig. 
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0.3 



0.6 



m. 



Normalized differential decay probability for r] ir^jj as a function of 
/m^. The dotted line represents the phase space distribution. The dashed 
line corresponds to the calculation at 0{p^) without reduction of the resonances. The 
dash-dotted lines display the uncertainty in the reduction of meson resonances for differ- 
ent values of the parameter Z\. The solid line corresponds to the values of the coefficients 



Li and di used in Ref. 



77 7r°7r° cross section results from the interference between the Born amphtude 
on the one hand and one- and two-loop amplitudes on the other hand. Thus the ex- 
perimental data are not sensitive enough to the various Born contributions described 
by di. On the other hand, the Born contribution is dominating in the r] n'^^f^f 
decay. In Fig. 2 we show the normalized differential decay probability as a function 
of m^^/m^. In this case the differential distribution is very sensitive to the input 
parameters di. Thus data of the differential distribution would be of great value for 
constraining these parameters. 

In conclusion, a self-consistent, quantitative description of 77 — > n^n^ and f] — > 
7r'^77 data at 0{p^) is still problematic. A good description of the 77 tt^tt^ cross 
section has been achieved whereas a satisfactory, quantitative prediction of the decay 
width seems to be beyond the reach of an ordinary calculation at 0{p^). 
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